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During the last six months, we have investigated the transverse heating of 
H+ ions and the consequent dynamics of the polar wind. This work is reported 
in a recent paper [Singh, 1992]. A copy of the paper is enclosed. We have also 
examined the effect of transverse ion heating using a particle code. This work is 
in progress. Some initial results were presented in two papers at the last IUGG 
meeting in Vienna [Singh, 1991 a,b]. 
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plasma perturbations Created by transverse ion heating 
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Transversely healed ions provide an important mechanism for the transport of ionospheric ions to the 
distant parts of the magnetosphere. The heating occurs via wave-particle interactions. The transport 
of such ions has been generally treated by a test particle approach, according to which ions move upward 
adiabatically under the action of the mirror force. The perturbations in the plasma created by the 
heating are studied here for the first time by solving a self-consistent set of plasma transport equations, 
in which heating is included phenomenologically. The nature of the perturbation is found to depend on 
the field-aligned extent of the heating region and on the heating time. When the heating occurs in a 
relatively localized region, the perturbation consists of a density cavity topped by a density enhancement 
(bump). The plasma in the bump is that expelled from the cavity. For weak heating the bump 
propagates upward without any steepening near the density fronts. For relatively strong heating the 
bump steepens into a forward-reverse shock pair. The forward shock (F) forms near the leading edge 
while the reverse shock (R) is at the top of the density cavity, near the trailing edge of the bump. When 
the R-F shock pair has formed, the velocity profile shows double "sawtooth" structure. Such shock 
structures are similar to the features of interplanetary shocks associated with solar wind transients. 
Parallel temperatures of the ions are reduced inside the cavity and increased in the density bump 
adiabatically. The transversely heated ions occur inside the upper part of the cavity and in the density 
enhancement. The dynamics of the transversely heated ions is not determined by the mirror force alone. 
The electric fields developing in the perturbation significantly affect their upward motion; the electric 
fields near the lower edge (R) of the density bump are downward, and they retard the upward motion. 
The electric and mirror forces on the transversely heated ions are compared. The electric fields near 
the leading edge (F) of the density bump are upward, and when they are sufficiently strong they 
compress the ions in the ambient polar wind. This creates an ion population heated only in T f just near 
the leading edge of the perturbation. During the course of its evolution the electric field distribution 
shows pulselike features. For impulsive heating the distribution evolves into a predominantly unipolar 
upward pointing electric field near the leading edge. Such an electric field pulse moves upward and its 
field strength depends on the electron temperature and the density gradient across F. The perturbation 
produced by transverse heating extending over long distances is characterized by an upward expanding 
density cavity with its lower end anchored at the lower edge of the heating region. Inside the cavity the 
spatial distributions of density, flow velocity, and parallel and perpendicular temperatures reach a 
quasi-steady state while above it they continue to evolve. The dynamics of the ions in the cavity, 
including their pitch angle distribution, is described by a test particle approach under the influence of 
the mirror force. When the heating ceases, the cavity fills from the bottom through the expanding polar 
wind, and the transversely heated ions inside the cavity ride on top of this polar wind expansion. 


1. INTRODUCTION 

During the last decade, transversely heated ions have been 
extensively studied [e.g., Klumpar, 1986], The transverse 
heating occurs by wave-particle interactions (see, for example, 
review by Ashour-Abdalla and Okuda [1983]). Both localized 
heating and extended heating have been studied in the 
literature. In the case of localized heating, the ions are heated 
in a region of limited extent, in which ions attain a pitch angle 
near 90°. As the ions move upward under the influence of the 
mirror force in the northern hemisphere, their pitch angle 
increases adiabatically (e.g., Klumpar , 1979], When the heating 
is extended [Temerin, 1986], ion continually gain energy, and 
their pitch angle attains an asymptotic value [ Change t al., 1986; 
Singh and Hwang, 1987], The upward motion of the heated 
ions is generally considered by use of a simple test particle 
approach. The perturbations in the plasma created by the 
transverse heating, including its effect on the distributions of the 
plasma density, ion temperatures, and flow velocity, have not 


been studied so far by solving a self-consistent set of plasma 
equations. 

The purpose of this paper is to report our first attempt to 
understand the dynamics of the plasma perturbations created 
by transverse ion heating by solving an appropriate set of 
plasma transport equations similar to those for the polar wind. 
It is found that the heating-induced perturbations in the density, 
velocity, and temperature profiles of the ambient polar wind are 
quite severe. The temporal evolution of perturbations created 
by localized and by extended heating is presented here. For 
localized heating, the perturbation consists of a plasma cavity 
topped by a density bump which eventually evolves into a 
reverse-forward shock pair if the heating rate is sufficiently 
large. The parallel temperature structure in the perturbation 
is found to depend on the magnitude of the heat flow. Heat 
flow in a collisionless plasma is not a well-understood 
phenomenon. Therefore in the hydrodynamic model adopted 
in this paper, heat flows are included phenomenologically. 
Anomalous plasma effects may severely inhibit the heat flux 
[Manheimer, 1977]. The parallel temperature distributions for 
different heat flow conditions are presented. 

The density cavity and the bump are the results of the 
expulsion of the heated ions by the mirror force. Thus it is 
expected that the ions in the entire bump be transversely 
heated. This is found to be true only initially. As the bump 
evolves, its leading edge compresses the ambient polar wind as 
a snowplow compresses snow ahead of it. This newly 
compressed plasma becomes part of the bump, but it is heated 
only in Ty. Consequently, the transversely heated ions are 
found near the top of the cavity and in the lower part of the 
density bump, and not near its leading edge. However, the 
exact topology of the ions with enhanced Ty and/or T ± 
depends on several factors, such as the heating time and the 
electron temperature. 

The electric field distributions associated with the 
perturbations in the plasma flow show pulselike features. The 
electric fields in the pulses can point downward or upward 
depending on their location with respect to the density bump. 
When heating is very brief in time and localized in space, the 
evolving perturbation has an electric field distribution 
dominated by an upward pointing electric field pulse near its 
leading edge. Such fields appear as upward moving weak 
double layers. 

Electric fields developing inside the perturbations are found 
to have significant effects on the upward motion of the 
transversely heated ions. Just before the reverse shock forms, 
the downward electric field force near the trailing edge of the 
bump exceeds the mirror force on the ions. At later times the 
two forces are found to be comparable. Thus the upward 
motion of the transversely heated ions is determined not by the 
mirror force alone, but by the self-consistent evolution of the 
entire perturbation. 

When the heating is extended, the major feature of the 
perturbation is the formation of a deep plasma cavity topped by 
an extended density enhancement, which is relatively weak. The 
cavity progressively expands upward. Inside the density cavity 
the flow reaches a quasi-steady state which is maintained as 
long as the heating lasts. The relaxation of the polar wind after 
the heating ceases is presented. 


Previously, time-dependent polar wind studies have been 
carried out by several authors: Singh and Schunk [1985] 

studied the evolution of reverse and forward shock pairs when 
the polar wind was artificially disturbed by density cavities and 
bumps. The present study shows how such cavities and bumps 
are naturally created by transverse ion heating. Gombosi and 
Killeen [1987] studied transients in response to low-altitude ion 
heating. Mitchell and Palmadesso [1983] and Ganguli and 
Palmadesso [1987] have studied the time-dependent features of 
the polar wind in response to the wave-particle interactions 
driven by field-aligned currents. The latter authors have dealt 
with the issue of transverse ion heating, but the features of the 
perturbations in response to transverse heating of the polar 
wind, as seen in the present paper, have not been reported. 

The rest of the paper is organized as follows; in section 2 the 
theoretical model is described. The results of the perturbation 
and its evolution are described in section 3, which includes 
discussions of the effects of heating rate, heating time interval, 
heating region size, electron temperature, and heat flow on the 
evolution. The paper is concluded in section 4. 


2. THEORETICAL MODEL 


We model here the transport of plasma in response to the 
transverse heating of the polar wind at altitudes of >3000 km. 
The relevant equations governing the magnetic-field-aligned 
transport of the mass, momentum, and thermal energies of 
hydrogen ions in a collisionless plasma are given by 
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where t is time; r is the geocentric distance along the flux tube; 
s is the distance along the tube from the lower boundary of the 
flux tube at r = r Q (see Figure 1); V t , T^, and Tj are the 
number density, flow velocity, and parallel and perpendicular 
temperatures of the polar wind, respectively; E is the parallel 
electric field; g| is the component of the gravitational force 
parallel to the magnetic field; and m,- and q i are the ion mass 
and charge, respectively. A is the cross-sectional area of a flux 
tube, and it is assumed that Aar 3 . Except for the time- 
independent terms, equations (1) to (4) are similar to the 
formulation given by Holzer et ai. [1971], who derived these 


equations directly from the moments of the Boltzmann 
equation. We do not solve the heat flow equations, which have 
proven to be quite troublesome to solve numerically 
[Palmadesso et al., 1988]. Instead, we have included the effects 
of heat flow heuristically by closely following the treatments in 
solar wind studies [e.g., Metzler et al., 1979]. 

In a collisionless plasma the usual picture of heat flow, given 
by q a = -K a V T a with K a as the thermal conductivity, may not 
be valid because X, the mean free path, is > >L T = (7~* 37V 
dr)' 1 , the scale length in the temperature variation. In such a 
situation, maximum heat flux is given by the transport of the 
thermal energy density ( nk B T a ) by the parallel thermal velocity 
Vth ~ where T a denotes T ± or Ty However, 

plasma turbulence has been found to inhibit the heat flow 
drastically [Manheimer, 1977]. Accordingly, we assume that 

<lcu = € 1 ' n i kT <xi v th ( 5 ) 

where c = -1 if 3TJ& > 0 and e = 1 if 3TJds < 0. Thus 
temperature gradients determine only the sign of the heat flow 
and not its magnitude. The factor tj is the reduction due to the 
anomalous plasma effects. In the present macroscopic 
calculations we cannot determine the value of r> self- 
consistently. Therefore, in order to develop a feel for the 
effects of heat flow, calculations are performed by varying i) 
from 0.01 to 1. 

The electric field E is obtained from the electron momentum 
equation by making the following assumptions: (1) plasma 
remains quasi-neutral, n e - (2) electron temperature (T e (s)) 

is isotropic and does not vary with time during the evolution of 
the plasma perturbation; and (3) electron flow velocity is 
calculated from the condition of no field-aligned current, i.e., 
n e V ( = n i V i or V g = Vj when n i = n e . Thus the electric field 
is given by 
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In the calculations presented here, we have limited the 
perpendicular ion heating to a relatively low level, for which it 
is found that the first two terms in equation (6) make negligibly 
small contributions to E. Since we do not solve the electron 
energy equation, care has been exercised to ensure that the 
spatial variation in the flow velocity has a negligible effect on 
the electron temperature. This is possible when the electron 
thermal velocity is considerably greater than the flow velocity. 
Calculations are presented with both constant and variable 
electron temperatures with altitude. 

It is likely that the electrons may develop a temperature 
anisotropy [Mitchell and Palmadesso, 1983]; its contribution to 
the electric field is proportional to 3 It, where r is the geocentric 
distance; inside the perturbation where the density and/or 
velocity have relatively large gradients, this contribution by the 
electron temperature anisotropy is found to be negligibly small 
compared to the third term in equation (6), unless the 
perpendicular electron temperature becomes comparable to the 
perpendicular temperature of the heated ions. Further 
quantitative comparison between the various contributions to E 
will be given later in the paper. 


Equations (1) to (4) are written in the form of a generalized 
continuity equation and are solved numerically by the flux- 
corrected transport technique [Boris and Book, 1976]. These 
equations are solved as an initial boundary value problem. The 
initial conditions correspond to the steady state (ambient) polar 
wind, in which protons flow supersonically. The boundary 
conditions are as follows; at the lower boundary (r = r Q ) we 
assume that n = N 0 , V - V Q , and T r = T\ = T 0 , where N a , 
V , and T 0 are the ionospheric density, flow velocity, and 
temperature at the altitude r = r Q . In the present calculation 
it is assumed that N a - 500 cm' 3 , T Q = 0.3 eV, and the value 
of V Q is varied from case to case depending on the electron 
temperature T e , which exercises a great deal of control on the 
flow velocity. At the top boundary (r = r max ) the continuity of 
the flow is maintained. 

The heating rate STJSt appearing in (4) is given by 

(ST ± /St) lmb - 0.25 ( 7 ) 

where i|i is the wave power spectral density in units of V 2 m' 2 
Hz' 1 . Such a heating rate is the consequence of the interaction 
of ions with broadband waves when the wavelengths are much 
larger than the ion Larmor radius [Singh and Schunk, 1984; 
Chang et al., 1986]. It is worth mentioning that knowledge of 
the exact nature of the heating is not important here, as we are 
studying the large-scale effect of the heating on the polar wind 
rather than the microscopic causes of the heating. In the 
calculations presented here, different values of t|r are adopted 
depending on the desired heating level. 

3. PLASMA PERTURBATION 

We present here the nature of the plasma perturbation 
created by the transverse heating of the ions. The temporal 
evolution of the spatial structure of perturbations for different 
heating rates, size of the heating region, and electron 
temperature distributions is discussed. 

3.1. Localized Heating 

It is instructive to examine the nature of the perturbations for 
a relatively low heating rate, which produces relatively weak 
perturbations. Figures 2a to 2d show the evolution of the 
perturbation when the heating occurs over a distance As h = 210 
km starting at a height of 5100 km (s = 2100 km). The 
heating rate is 0.24 eV/s, which is obtained from (7) for a 
power spectral density i|r = 10~ 8 V 2 m~ 2 Hz -1 . The heating 
time interval A/ h is 1.844 min. The electron temperature 
distribution is assumed to be a constant, T e (s) = 1 eV. The 
heat flow reduction parameter v is 0.5. Figures 2aio2d show 
the evolution of density, flow velocity, parallel temperature, and 
perpendicular temperature, respectively. The profiles shown in 
this figure are at times t 0 , t v t 2 , r 4 , and t 6 where r n = n x 0.922 
min. The profiles at t = t 0 correspond to the initial 
unperturbed polar wind. The perturbations in n, V, and T| are 
minimal at = 0.922 min as they are hardly visible at the 
scales used in the figures. However, at this time the 
perturbation in T indicates a heating up to a maximum value 
of -3 eV. Eventually, the perpendicular heating affects other 
quantities, and the perturbations in n, V, and Tj can be clearly 
seen for t > / 2 . 


The perturbation in the density consists of a cavity topped by 
a density enhancement. This perturbation is simply the 
consequence of the expulsion of the plasma by the mirror force 
acting in the region with enhanced T L . The density 
perturbation is accompanied by an enhancement in the flow 
velocity (Figure 2b) and by a reduction in the parallel 
temperature inside the plasma cavity and an increase inside the 
density enhancement. The structure in the parallel temperature 
is accounted for by an adiabatic expansion inside the cavity and 
a compression inside the density enhancement [e.g., 
Tanenbaum, 1967]. For example, the minimum parallel 
temperature (T| min ) in the plasma cavity and the maximum 
temperature (T. max ) in the density enhancement are found to 
approximately obey the relation 

^'|min^ r |max = (^min^max) ^ 

which comes from the adiabatic law = const with 

Y = 3 for a one-dimensional process. 

Figures 2a to 2d show that as long as the heating continues 
the lower end of the perturbation remains fixed at the altitude 
of the lower edge of the heating region. This is seen very 
clearly by comparing the profiles at t — tj and f 2 in Figure 2d. 
When the heating stops, the entire perturbation moves upward. 
The upward moving perturbation rides on the upward flowing 
polar wind. The peak in the density enhancement is seen to 
move with a velocity V n = 14 km/s, which is considerably 
slower than the flow velocities inside the perturbation. As the 
transversely heated ions move upward after the heating stops, 
the maximum value of T x continually decreases. This decrease 
is found to be approximately adiabatic, i.e., 7, maT a B. 

EnJianced heating rate. When the heating rate is enhanced, 
the evolution of the plasma perturbation is seen to change 
considerably. Figures 3a to 3d show the evolution of the 
perturbation when the heating rate ( STJdt ) is 2.4 eV/s, in 
contrast to 0.24 eV/s for the perturbation shown in Figures 2a 
to 2d. All other parameters for Figures 3a to 3d are the same 
as for Figures 2 a to 2d. When t = fj = 0.922 min, the 
transverse heating has increased 7" to a maximum value 
T max = 30 eV (Figure 3d), which is about 10 times the 
maximum T ± at / = tj in Figure 2d. As in Figures 2 a to 2c, 
Figures 3a to 3c show the formation of a plasma cavity topped 
by a density bump, enhancement in the flow velocity throughout 
the plasma perturbation, a decrease in the parallel temperature 
in the cavity, and an increase in the bump. However, there are 
significant differences between the magnitudes of the 
perturbations shown in Figures 2a to 2d and Figures 3a to 3d. 
Consequently they evolve quite differently. 

The large-amplitude perturbations in density and velocity 
profiles produce nonlinear effects. A prominent effect of 
nonlinearity is that the leading edges of the perturbations in 
Figures 3a to 3c are relatively much steeper than those seen in 
Figures 2 a to 2c. This is a common feature of the nonlinear 
evolution of ion acoustic shocks and solitary waves [Tidman and 
Krall, 1971]. 

The leading edge of the perturbation indicated by F in 
Figures 3a to 3c is a forward shock, which connects the 
perturbation with the ambient polar wind. The forward shock 
(Kfs) velocity is calculated by the conservation of mass flow by 


the relation 


^fs = 


"1 V \ ~ n a V a 


(9) 


where and V 1 are the peak density and velocity, respectively, 
near the forward shock F and n a and V a are the corresponding 
quantities of the ambient polar wind just above the shock. At 
t = t 6 , = 0.178 n Q , K x = 48 km/s, n a = 0.047 n Q , and 

K a = 29 km/s, giving ^fs = 32 km/s, which is found to be in 
agreement with the shock velocity calculated from the time of 
flight; the shock positions at t - t 5 and t 6 are at s - 7560 km 
and s - 9240 km, respectively, giving Vfs -31 km/s. This 
agreement indicates that the shock is quite steady as it moves 
upward. 

There is an interesting process which occurs near the leading 
edge (F). The perturbation in parallel temperature runs slightly 
ahead of that in the perpendicular temperature. This is shown 
by the dashed vertical line running through all four panels of 
Figure 3. On the right-hand side of this line the ions in the 
density enhancement are not heated in T ; instead, they are 
heated in T. up to about 2 eV. It is worthwhile to understand 
the origin of such ions in the perturbation. Initially, when the 
transverse heating produces a perturbation in the form of a 
density enhancement by expelling plasma via the mirror force, 
the entire plasma in the perturbation is heated in T . For 
relatively strong heating, as shown in Figure 3, the density 
enhancement is strong. The upward electric field and the 
pressure forces near the leading edge (F) are able to further 
compress the plasma in the ambient polar wind just ahead of 
the perturbation produced by the mirror force. The newly 
compressed plasma is not heated in T v but the adiabatic 
compression leads to the enhancement in Ty As will be shown 
later, this effect is considerably enhanced when an increase in 
the electron temperature enhances the electric field near F (see 
Figure 10). This process of ion heating has been seen 
previously in plasma expansion studies [Mason, 1971; Singh and 
Schunk , 1983], and more recently it has been termed an 
electrostatic "snowplow" [Schunk and Szuszczewicz, 1988]. 

Just like the pressure and electric forces near the leading 
edge, the downward forces near the trailing edge of the density 
enhancement significantly affect the evolution of the 
perturbation. In particular, such forces lead to the formation 
of a reverse shock (R) near the trailing edge of the density 
bump. At an early time ( t < f 2 ) the velocity profile (Figure 36) 
has only one "sawtooth" distribution. At later times (/ > t 4 ) it 
develops two sawtooth distributions, one near the leading edge 
(F) and the other near the trailing edge. The latter sawtooth 
distribution begins to form at about t - t 4 . This sawtooth 
formation proceeds as follows. When the downward electric 
and pressure forces associated with the sharp gradients near the 
trailing edge of the density bump become strong enough, the 
flow is slowed down there. Such local slowing down of the flow 
leads to the formation of the new sawtooth in the flow velocity. 

The structure in the flow near the new sawtooth is a reverse 
shock (R), where the upward jumps in the density and T| are 
accompanied by a downward jump in the velocity. In contrast 
to this both velocity and density jumps are downward across the 
forward shock. 


Reverse shocks are known to form in a variety of situations, 
such as the transient solar wind [Hundhausen, 1973] and some 
types of nozzle flows [Liepmann and Roshko, 1957], The 
physical factor facilitating the formation of a reverse shock is 
the creation of a low-pressure zone in a high-speed flow. 
Transverse heating produces such a situation when the plasma 
is evacuated by the mirror force. However, it is important to 
note that only when the heating and the associated 
perturbations in the plasma are strong enough, do the reverse 
and forward shocks form. This point is well illustrated by 
comparing the results in Figures 2a to 2d with those in Figures 
3 a to 3d, which are for weak and strong heating, respectively. 

The evolution of the electric fields during the evolution of the 
plasma perturbation is shown in Figure 4. The electric fields 
are primarily determined by the slopes in the density (dn/ds 
term in equation (6)) profiles. The time variation of the 
velocity and the slopes in its profile do not significantly 
contribute to the fields. The solid curve in Figure 4 is the 
contribution by dn/ds alone while the circles show the sum of 
the contributions of all three terms in equation (6) for the 
electric field. It is clearly seen that the two curves are 
remarkably close, implying that the contributions from the 
terms containing dV/dt and dV/ds in equation (6) are 
insignificant. 

The role of the downward electric field in the formation of 
the reverse shock can be appreciated by noting that at t = / 4 a 
large downward (negative) electric field develops. The 
maximum value of the negative electric field is about -9 <iV/m. 
This field plays a crucial role in slowing down the flow and 
thereby in the formation of the sawtooth "R," where the reverse 
shock forms. In order to clearly demonstrate this point, we 
have compared in Figure 5 the mirror force on the ions with 
the electric force arising from the electric fields associated with 
the perturbations. The mirror force F m is calculated using 

3 T ■ 

F = e (10) 

m r 

where T Ll is in electron volts and r is the geocentric distance. 
The solid and dashed curves in Figure 5 show the electric field 
and F m /e in units of microvolts per meter, respectively. We 
note that at t s f 2 > F m exceeds the downward electric force eE 
almost over the entire perturbation. However, at t = r 4 the 
maximum downward electric field is 9 <iV/m while at the same 
location the upward mirror force per unit charge is <4 iiV/m. 
The excess of the downward electric force locally slows the 
transversely heated ions forming the sawtooth in V at R. At 
later times this trend is continued as T decreases with 
increasing altitude. 

Effects of heating duration. It is certain that the transverse 
heating events in space occur over a variety of temporal and 
spatial scales. Here we discuss the nature of the perturbation 
when the heating time interval (A/ h ) is reduced while the 
heating rate is enhanced to achieve a desired heating level. 
Figure 6 shows results for Ar h = 0.1844 min and (dT^/ar),^ 
- 24 eV/s, while As h = 210 km and r? = 0.5 as in Figures 2 and 
3. The basic nature of the perturbation in this figure is the 
same as in Figure 3. However, there is a difference in its 
evolution. We noted earlier while discussing Figures 2 and 3 


that the lower end of the perturbations remains fixed as long as 
the heating is on. When the heating is over, the composite 
perturbation moves upward, riding on top of the ambient polar 
wind expanding upward. Since the heating time is 0. 1844 min 
in Figure 6 in contrast to 1.844 min in Figure 3, the composite 
perturbation begins to ride on top of the upward expanding 
polar wind earlier in Figure 6 than in Figure 3. Consequently, 
the shorter the heating event, the higher the perturbation 
moves upward. Comparing Figures 3 a and 6a, we find that in 
the latter figure the leading edge of the perturbation has 
reached the distance s = 12,000 km in contrast to 10,000 km 
in Figure 3. 

Effect of heat flow on the perturbation. In the above 
discussion we assumed that the heat flow reduction factor tj is 
0.5 in equation (5). Since the value of t? is uncertain, it is 
worthwhile to carry out a parametric study that varies t> to see 
how the variation affects the evolution of the plasma 
perturbation. Figure 7 shows the comparison of the 
perturbations for i? = 0.01 (dashed curves) and jj = 1 (solid 
curves), respectively. The panels in the left- and right-hand 
columns show the perturbation at early and late times, 
respectively. The choice of ij = 0.01 implies a negligibly small 
heat flow. On the other hand, »j = 1 yields the maximum heat 
flow. The early time perturbation in the left-hand column 
clearly demonstrates that the heat flow has negligibly small 
effect on the density and flow velocity while temperature 
profiles are slightly modified. On the other hand, at relatively 
late times when the reverse shock has formed, even the density 
and flow velocity profiles show significant modification by the 
enhanced heat flow. The most significant effects in the n(s) and 
V(s) profiles are found near the reverse shock, where flow 
velocity is relatively reduced. With reduced flow velocity the 
transport of the thermal energy due to the thermal motion 
becomes more important. The enhanced heat flow is seen to 
stretch out the region of enhanced parallel and perpendicular 
temperature. The density depletion and the reduction in the 
parallel temperatures in the plasma cavity are less severe with 
the enhanced heat flow. 

The heat flow also has significant effects on the propagation 
of the reverse and forward shock pair. The forward shock is 
seen to move faster with increasing heat flow. On the other 
hand, the reverse shock slows down. Consequently, the density 
enhancement and the associated features in Tj and velocity 
distributions are stretched out. However, it must be said that 
the basic features of the perturbations created by the transverse 
ion heating are not severely affected by the heat flow. This is 
true at least for the time span over which the evolution for ij = 
0.01 and 1 is shown. 

Effect of increasing electron temperature with altitude. In the 
calculations presented so far, electron temperature T e was 
assumed to be constant all along the flux tube. Since electrons 
are expected to be progressively warmer with increasing 
altitude, it is useful to examine the effects of such electron 
temperature distributions on the evolution of the plasma 
perturbations. Not knowing an exact distribution, we have 
arbitrarily chosen here that the electron temperature linearly 
increases from 0.3 eV at s = 0 to 5 eV at s = 16,500 km. The 
choice of the temperature at the high altitude is guided by 
observations that at altitudes of about 1 Rg the thermal 


electron temperature can be up to several electron volts [e.g., 
Kintner et al., 1979]. 

Figures 8 and 9 show the evolution of the plasma 
perturbations for this electron temperature distribution. Other 
parameters for the calculations shown in these figures are the 
same as those for Figures 3 and 4. The basic features of the 
perturbation, namely, the formation of a plasma cavity topped 
by a density enhancement and its eventual evolution into a 
reverse-forward shock pair, are similar to those in Figures 3 a to 
3d. So is the evolution of electric field perturbations for the 
two sets of calculations with a constant electron temperature of 
1 eV (Figure 4) and with the electron temperature increasing 
with altitude (Figure 9). However, there are some interesting 
quantitative differences. For the assumed electron temperature 
distribution in Figure 8, T e (s) > 1-eV for the entire 
perturbation in contrast to a constant 1 eV temperature over 
the entire perturbation shown in Figure 3. Thus, wherever a 
sharp density gradient develops in the perturbation, the 
magnitude of the electric field becomes larger for the 
perturbation in Figure 8 than for the perturbation in Figure 3. 
For example, at t = r 2 , the electric field (Figure 9b) near the 
trailing edge of the density bump (Figure 8a) is nearly twice the 
corresponding electric field (Figure 4b) with the constant 1-eV 
value of T e . This hastens the formation of the reverse shock 
near the trailing edge in Figure 8. Since the reverse shock 
moves slower than the forward shock, the density bump (Figure 
8a) and the associated enhancement in the parallel temperature 
(Figure 8c) are considerably stretched out. This stretching of 
the perturbation can be clearly seen by comparing the electric 
field distributions at t - t 6 in Figures 4 d and 9 d. 

The electric field enhancement near relatively large density 
gradients caused by an increase in the electron temperature is 
found to affect the upward motion of the transversely heated 
ions. For example, Figure 8 d shows that at t = t 6 = 5.6 min, 
transversely heated ions have traveled upward to a height s - 
8200 km in contrast to a considerably greater height of 9300 
km in Figure 3d, in which the electron temperature has a 
constant value of 1 eV. The heights quoted here are for the 
locations at which TJs) suddenly begins to fall with increasing 
altitude. This retardation of the transversely heated ions is a 
direct consequence of the enhancement of electric fields due to 
enhanced electron temperature. As mentioned earlier, the 
downward electric fields retard the motion of the upward 
moving ions. The downward fields in Figure 9 are found to be 
consistently higher than those in Figure 3. 

Perturbations associated with impulsive ion heating. In the 
previous examples of ion heating, the heating was relatively 
extended over a distance of As h = 210 km. We now wish to 
examine the nature of the perturbation when the size of the 
heating region is reduced. In the example presented here we 
assume As h - 60 km. Owing to numerical constraints we are 
not able to further reduce As h and at the same time follow an 
extended length of the flux tube to allow for the evolution of 
the perturbation. Along with the size As h , the heating time Ar h 
is also reduced to 5 s. Thus the heating simulated in this 
example is quite impulsive in both time and space. In order to 
have a substantial perpendicular heating for the limited values 
of As h and A/ h , the heating rate is enhanced to 240 eV/s, which 
corresponds to iji = 10' 5 V 2 m' 2 Hz' 1 in equation (7). We 


have assumed here that electron temperature T e is 10 eV over 
the entire region of the perturbation and that tj = 0.5. 

Figures 10 a to lOd show the evolution of the perturbations up 
to t = t 4 = 2 min. Note that t n = n x 30 s in this figure. As 
in previous examples, we find that at an early time (t < tj) the 
basic feature of the perturbation is the formation of a plasma 
cavity topped by a density enhancement. However, at later 
times the evolution differs from the previous example in a 
significant way concerning the formation of the reverse shock. 
The density and velocity profiles show that the reverse shock 
does not develop fully. The sawtooth in the velocity profile 
near the trailing edge of the bump remains weak all along the 
evolution shown in Figure 106. It is seen that the leading edge 
(F) of the density bump is quite steep, but its trailing edge is 
seen to smooth out with increasing time. This is attributed to 
the fact that in the impulsive heating lasting over a short 
duration, the cavity begins to refill from the bottom at an early 
stage through the upward expansion of the polar wind. Thus 
the pressure inside the cavity does not become low enough for 
the reverse shock to develop fully. The forward shock and its 
dynamics are found to be quite similar to those in the previous 
case as discussed in connection with Figure 3. 

The electric field perturbations for the impulsive heating are 
shown in Figures lOe to 106. The most noteworthy feature of 
the electric field distribution shown in this figure is its evolution 
to a predominantly unipolar upward pointing electric field near 
the leading edge of the density bump when t>t 2 - 1 min. This 
dominant electric field pulse propagates upward with a velocity 
of about 60 km/s, which is about twice the H + ion acoustic 
speed with 10-eV electrons. Such upward propagating electric 
field pulses appear quite similar to the predominantly unipolar 
electric fields observed in the auroral plasma [ Temerin et at., 
1982]. However, the observed fields are generally interpreted 
as ion acoustic double layers which have scale lengths of a few 
tens of Debye lengths (<100 m). In contrast, in the present 
calculations we have a spatial resolution of -60 km, which is the 
minimum intergrid spacing used in our calculations. Therefore, 
if there is any connection between the observed electric field 
pulses and those pulses generated by impulsive ion heating 
presented in this paper, it remains to be explored. 
Observationally, this task requires a close scrutiny of the plasma 
data along with the data on electric fields, correlating the 
perturbation in the plasma density with the electric field pulses. 
Theoretically, it is essential that the temporal and spatial scales 
of the calculations be reduced to the extent that the electronic 
field pulses of durations of a few milliseconds be produced. As 
mentioned earlier, such temporal resolution becomes extremely 
demanding in terms of computer resources. Therefore the 
main point of the discussion here is to suggest that 
perpendicular ion heating can produce parallel electric fields. 
We cannot rigorously compare such fields with double layers. 

Figures 10a to 10c again show that the leading edge of the 
perturbation travels ahead of the transversely heated ions. As 
discussed earlier, the primary reason for this behavior is that the 
large electric field near the leading edge acts like an 
electrostatic "snowplow" compressing the plasma just ahead of 
it. This compressed plasma ahead of the transversely heated 
ions has enhanced Tj but not T ± . This is clearly seen from 
Figures 10c and lOd. As in Figure 3, the dashed vertical line 


separates this compressed plasma from the transversely heated 
ions. 

3.2. Extended Heating Events 

A role of extended transverse heating of ions in producing ion 
conics having nearly constant cone angles with variation in the 
geomagnetic latitude has been suggested in explaining 
observations from the ISEE 1 satellite of low-energy ion conics 
in the outer region of the plasmasphere [Horwitz et al., 1982; 
Singh, unpublished manuscript, 1980]. More recently, extended 
ion heating has been proposed as a mechanism for producing 
ion conics in the auroral plasma [Temerin, 1986], in the plasma 
sheet [Chang et al., 1986], and in the refilling zone of the outer 
plasmasphere [Singh and Hwang, 1987]. However, all these 
studies have been essentially based on a test particle approach, 
in which the plasma perturbation is ignored. We discuss here 
the nature of the plasma perturbation produced by transverse 
ion heating occurring over an extended region along an open 
flux tube. Figures 11a to lid show the evolution of the 
perturbations in density, velocity, parallel temperature, and 
perpendicular temperature profiles, respectively, in response to 
an extended heating event. The perturbations shown in this 
figure are for a relatively low heating rate of 0.24 eV/s above 
the distance s = s Q = 6000 km (altitude 9000 km). The 
electron temperature is assumed to be 1 eV all along the flux 
tube. We assumed that tj = 0.5. Qualitatively similar results 
were obtained when T e was allowed to vary with altitude 
according to T e (s) = 0.3 eV + 4.3 x 10" 5 s eV/km, which yields 
a temperature of 0.3 eV at the boundary s = 0 and 2.5 eV at 
the top of the flux tube considered in the calculations (Figure 
l). In Figure 11, f n = n x 1.844 min. The evolution of the 
density profile (Figure 1 la) shows that the primary effect of the 
heating is the evacuation of plasma beginning at s = s Q . The 
plasma cavity formed by the evacuation continually expands 
upward. The expansion velocity of the cavity front after the 
initial stage (/ * r 3 ) is found to reach an asymptotic value of 50 
km/s. 

Figure 11a shows that density above the expanding cavity 
increases only slightly in comparison to the increase in the 
density enhancements for localized heating. This is primarily 
due to the continual upward acceleration of the ions by the 
mirror force resulting from continual transverse heating in 
space. Acceleration of ions resulting in a velocity increment 
above the cavity can be seen from Figure 1 lb. It is found that 
at any given time the velocity and density profiles reach 
quasi-steady states below the leading edge of the cavity. Above 
the leading edge of the density cavity the flow velocity is still 
evolving through the acceleration of ions by the mirror force. 
Inside the cavity, where the steady state is reached, the mirror 
force is found to balance the nonlinear dlds(l/2 V^) term in 
the momentum equation, i.e., 

4 (i ("“?>) -Wjr; 1 ") 

In order to demonstrate this force balance we compared the 
ion flow velocity obtained by integrating the above equation 
with the velocity profiles shown in Figure 116. T xi (r) used in 
the integration are the same as in Figure 11 d. At any time the 


integration is carried out from bottom to top of the cavity. The 
comparison is shown in Figure 12. The dashed curves in 
Figures 12a to 1 2d are the product of integrating (11), while 
solid curves are the same as in Figure 116. It is seen that 
equation (11) predicts the flow velocity inside the cavity quite 
accurately. In view of this close agreement, it is interesting to 
note that during the steady state phase the fluid flow is 
essentially given by the test particle approach. In this situation 
the pressure and electric forces are dominated by the mirror 
force. However, this is true only if the transverse heating 
results in relatively large perpendicular temperatures. 

The evolution of the parallel temperature (Figure 11c) shows 
that ions cool down inside the cavity while they are slightly 
heated above it as the plasma is compressed to a slightly higher 
density. Figure 1 Id shows that as the ions enter the heating 
region they quickly acquire a relatively large perpendicular 
temperature. This results in parallel acceleration of ions 
through the mirror force. Therefore, as ions move up, their 
transit time through a given distance in space continually 
decreases. This leads to a continually decreasing spatial 
increment in T , despite a constant heating rate with respect 
to time. This feature is borne out from Figure lid. The log 
scale used in Figure lid does not allow one to clearly see the 
evolution of T (s). Therefore T x (s) is replotted in Figure 13. 
It is seen that after the cavity forms at t > r 3 = 5.5 min, T is 
lower inside the cavity than above it. Inside the cavity, T ± (s) 
attains a quasi-steady profile as n(r), V(s), and T j(S) do. 

It is useful to examine the spatial evolution of the ion pitch 
angle defined by 

1/2 

/Vj ( 12 ) 

where V i is the parallel flow velocity. Figure 14 shows the 
distribution of a p (s) at / = / g , when the leading edge of the 
density cavity shown in Figure 1 la is just about to exit from the 
top of the calculation region which extends up to s = 52,500 
km. As expected, below the heating region, is found to 
steadily decrease with increasing height. Near tne lower edge 
of the heating region, a p abruptly attains a large value over a 
relatively short distance of about 2000 km. The maximum pitch 
angle is about 63° (or 1 17° in the northern hemisphere). Then 
it steadily decreases and asymptotically reaches 30°. This 
asymptotic value of the pitch angle can be analytically estimated 
as follows. Using a test particle approach, the motion of an ion 
through the heating region is described by [Singh and Hwang, 
1987] 

^ = k B {ZTJ3t ) turb Vj 1 (13) 
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where W = ZitnV\ + 2k B T x . When a p tends to its asymptotic 
value, datp/dr = 0, yielding 
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Assuming a p is nearly constant when daJdr -> 0, equation (13) 
is integrated to give W^cos a p = r(dT/dt) turb and 

tan a p - l/y/3" or a p - 30° ( 16 ) 

This is the asymptotic value of a p as seen from the pitch angle 
distribution in Figure 14. 

The above calculation shows that the behavior of the heated 
ions is quite well predicted by a test particle approach. But in 
the heating process the plasma is severely affected, and plasma 
perturbations cannot be described by this approach. 

The heating event shown in Figures 1 la to 1 Id lasts for 20 
min. After this period the plasma in the flux tube begins to 
recover as shown in Figure 15. The recovery is primarily 
through the upward expansion of the polar wind into the 
plasma cavity created by the transverse ion heating. The 
recovery seen here is somewhat similar to the polar wind 
expansion into a depleted flux tube previously studied by Singh 
and Schunk [1985] except for the presence of warm ions inside 
the cavity. Some interesting features of the recovery process 
are as follows. As the polar wind expands upward, its ambient 
properties are restored below a certain height. This height is 
seen to increase with a rate of about 8 km/s, which is 
comparable to the ion acoustic speed near the lower part of the 
flux tube. When t = r 24 = 44 min, the recovery has occurred 
only up to a height s = 15,000 km, as indicated by a big open 
square on the density profile curve in Figure 15. As the cavity 
refills from the bottom, the density in the upper part of the flux 
tube is seen to continually decrease. This is primarily due to 
the relatively large velocity flow in the upper part of the cavity. 
The transversely heated ions move upward on top of the 
expanding polar wind (Figure 15c/). Such ions in the plasma 
cavity are seen to progressively cool in T . The cooling is 
primarily due to heat convection. As T ± is lowered by cooling, 
the flow velocity tends to relax toward the ambient polar wind 
flow velocity (Figure 156). 

It is important to note that the maximum flow velocity 
attained in the extended ion heating event is about 140 km/s. 
This velocity is <20% of the electron thermal velocity along the 
flux tube. Thus the approximation made in calculating the 
electric field (equation (6)) is approximately satisfied. If the 
evolution with heating is to be followed until later times, the 
electric field expression needs further improvement by relaxing 
the isothermal assumption of the electron temperature. 

It is worth noting that the cavity on top of the expanding 
polar wind relaxes quite slowly. Following an extended heating 
event, the polar wind in the flux tube may take a few hours to 
relax. For example, in Figure 15 the ambient polar wind front 
reaches the distance r = 52,000 km (8 R^) in about 2 hours. 
During this recovery time the density profile in the flux tube 
shows two distinct regions; below the expansion front the 
density profile is that of the ambient polar wind while above it 
the density is generally lower than the ambient polar wind 
density. Satellite observations do indicate such transitions in the 
radial density profiles in the polar region of the magnetosphere 
[Persoon et al., 1983]. However, the reason for the density 


depletion at high altitudes remains unknown. It is likely that 
the extended transverse ion heating at high altitudes may create 
the transition. 

4. CONCLUSION AND DISCUSSION 

The main purpose of this paper is to study the nature of 
plasma perturbations created by transverse ion heating in the 
auroral and polar region of the magnetosphere. This study is 
motivated by the fact that transverse ion heating is a pervasive 
phenomenon. The creation of the perturbations by transverse 
ion heating and their evolution are studied by solving a set of 
plasma transport equations along open flux tubes for geocentric 
altitudes r > 9371 km/s. Both localized and extended heating 
events in space have been considered. Likewise, both short- 
and long-duration heating events have been studied. 

For localized heating, the main characteristics of the 
perturbation are as follows. The perturbation in the density 
consists of a density cavity topped by a density bump. During 
the early stage the plasma in the bump is that expelled from the 
cavity by the mirror force. The flow velocity is enhanced over 
the entire perturbation. The parallel temperature is lowered 
inside the cavity while inside the density bump it is enhanced. 
The reduction and increment in T, occur approximately 
adiabatically (T| a n 2 ). The transversely heated ions occupy 
the region near the top of the cavity and the bump. However, 
near the leading edge of the bump there are some ions which 
are heated only in T ( but not in T ± . This is attributed to the 
electrostatic snowplow effect. This effect is enhanced when the 
electric field near the leading edge is enhanced, for example, by 
enhanced electron temperature. As long as the heating lasts, 
the lower end of the perturbation remains anchored to the 
lower end of the heating region while the upper end of the 
perturbation moves upward. After the heating ends, the entire 
perturbation rides on top of the upward expanding polar wind. 

When the perturbations are relatively weak, produced by a 
relatively low transverse heating rate, the density and velocity 
profiles in the perturbations do not show any significant 
steepening. However, when the heating rate is sufficiently 
enhanced, the density and velocity profiles show the formation 
of a reverse-forward shock pair. The forward shock forms near 
the leading edge (top) of the density bump. It evolves through 
the steepening of a relatively large amplitude perturbation in 
the density. The reverse shock forms near the trailing edge of 
the density bump. 

The interesting feature of the electric field distribution 
associated with the plasma perturbation is that the fields appear 
as localized pulses, both pointing upward (near F) and 
downward (near R). These pulses propagate upward with the 
density perturbation. The magnitude of the fields increases 
with the electron temperature. The transversely heated ions 
occur near R where the downward electric field causes a 
downward force on them. Such electric forces are comparable 
to the mirror force on the ions. It is shown that the reverse 
shock forms when the downward electric force exceeds the 
mirror force forming a "sawtooth" in the velocity profile. 

When the heating is impulsive in both time and space, the 
electric field distribution eventually evolves into a unipolar 
upward pointing electric field near the leading edge. Such fields 


appear as weak double layers, such as are observed in the 
auroral plasma at altitudes near 1 [Tenierin et al., 1982]. 
However, it must be pointed out that if there is any connection 
between such observed double layers and the electric field 
produced by the impulsive ion heating, it remains to be 
explored. The observed electric fields have spatial scale lengths 
of a few tens of Debye lengths (-100 m). In contrast to this, in 
the calculations presented here we have resolved distances of 
the order of 60 km, which is the minimum intergrid sparing 
used in our calculations. For the purpose of a rigorous 
comparison, we need calculations which can resolve distances up 
to -10 m. On the other hand, observationally it is essential that 
electric fields be correlated with perturbations in the plasma 
density. Nevertheless, present calculations show that 
perpendicular heating generates parallel electric fields. It is 
worth noting that electric field measurements may be one way 
of verifying the presence of plasma perturbations created by 
transverse ion heating. 

The nature of the plasma perturbations produced by heating 
occurring over an extended region is found to be different from 
that produced by localized heating. The dominant feature of 
the density perturbation produced by extended ion heating is 
the formation of a cavity expanding ever upward as long as the 
heating lasts. Inside the cavity and above it the perpendicular 
temperature and the flow velocity are enhanced. As expected, 
the parallel temperature is reduced inside the cavity, and it is 
slightly enhanced above the cavity because of the enhanced 
density. Inside the cavity the profiles of density, velocity, and 
parallel and perpendicular temperatures reach a quasi-steady 
state and remain unchanged as long as the heating lasts. This 
steady state is achieved when the mirror force is balanced by 
the rate of charge of momentum due to the velocity gradients. 
Since the flows are highly supersonic, the thermal effects on the 
flow are found to be insignificant. Inside the cavity the ion 
motion is found to be governed by the test particle motion. 
The pitch angle distribution is found to reach an asymptotic 
value of 30° at large geocentric distances as theoretically 
predicted by a test particle approach [Singh and Hwang, 1987], 
This shows that the motion of transversely heated ions can be 
described by a test particle approach. But it is important to 
note that in the process of creating the transversely heated ions 
the plasma is strongly perturbed and the perturbations are 
beyond the scope of such an approach. 

The calculations presented here are based on the assumption 
that heating results in a bi-Maxwellian distribution and the 
entire ion population is affected. Such bulk heating has been 
observed from DE 1 [Moore et al., 1985], In situations where 
resonant ion heating affects only a fraction of ions, resulting in 
a complex velocity distribution function, the results presented 
can be of only limited validity. However, a rigorous treatment 
of resonant ion heating requires a kinetic model. It should be 
mentioned that when the transverse heating occurs over a 
broad range of energy, its evolution under the action of the 
mirror force gives rise to conics in the velocity distribution 
function. When such conical distributions evolve with increasing 
altitude, their parallel velocity spread increases. If this happens, 
the reduction of parallel temperature in the plasma cavity may 
be limited. This problem is difficult to treat using a fluid model 
for the plasma. 


However, a 16-moment description of the hydrodynamic 
equations [Barakat and Schunk, 1982] shows that the parallel 
temperature couples to the perpendicular temperature through 
a heat flow term given by 2 q x A' 1 34/ds [e.g., Ganguli and 
Palmadesso, 1987], The nature of localized heating indicates 
that T has relatively large gradients near both bottom and top 
edges of the region in which transversely heated ions exist (see 
Figure 3). The effect of heat flows across such gradients on the 
parallel temperature can be estimated by a dimensional analysis 
of the equation V { dT ^ds - 2q^A'^ M/ds, giving 
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where the heat flow parameters (u and t) are defined in 
connection with equation (5), V- is the flow velocity, F (h is the 
thermal velocity given by ( kTyjm |) 1/2 , A S is the spatial scale 
length over which T x varies, and R is the geocentric distance. 

For the sake of concreteness, let us discuss (17) in view of the 
heating event shown in Figure 3. Near the lower edge of the 
region of transversely heated ions the temperature gradient in 
T x is upward, giving a downward heat flow q x which implies 
e = -1. Thus the coupling of 7^ with T ± through heat flow 
produces an additional cooling of ions inside the plasma cavity. 
This cooling can be physically reconciled as follows. The ions 
carrying thermal (nkT x ) energy downward are acted on by the 
upward mirror force, and thereby the spread in their parallel 
velocity is reduced. 

Near the upper edge of the region of transversely heated ions 
the heat flow q x is upward, implying t - +1 and AT, > 0. 
Therefore the coupling between T. and T x near the upper 
edge yields an enhancement in T.. This agrees with the 
common sense that as the ions carry neat (nkT x ) upward, their 
parallel velocity spread increases through the action of the 
mirror force. In hydrodynamic models this increase in the 
velocity spread implies enhancement in the parallel temperature. 
The magnitude of this enhancement can be assessed from (17); 
at / = r 4 in Figure 3 the estimates for the various parameters 
appearing in this equation are V x « 50 km/s, - 10 km/s, As 
- 500 km, R - 10,000 km, and T x - 20 eV, giving AT ( - 1.2i> 
eV. 

Depending on the v (<1) factor, AT| due to the coupling 
between T| and T x through the heat flow q x is less than the 
enhancement in Ty due to the compression of the plasma in 
the density enhancement. Thus the coupling between T. and 
T x through perpendicular heat flow is not a dominating factor 
for the parameters of the heating events considered in this 
paper. However, it must be pointed out that in heating events 
producing much more energetic ions, ATj due to Ty-T x 
coupling through the perpendicular heat flow can be significant, 
and then it cannot be ignored in the model. In the present 
paper we have dealt with heating events producing T only up 
to a few tens of electron volts, for which it appears that our 
hydrodynamic model is adequate. 

The density depletion in the cavity created by transverse ion 
heating is found to be more than an order of magnitude below 
the density of the ambient polar wind. The minimum density 
in the cavity is found to be <1 cm" 3 . Persoon et al. [1988] have 
reported such low densities in the fine structures of the auroral 


plasma cavity. If numerous heating spots occur randomly 
distributed in time and over the auroral region, the combined 
effect of such heating will be the creation of an extended cavity 
like the auroral plasma cavity observed from satellites [Calvert, 
1981; Persoon et al., 1988]. The combined effect of numerous 
heating spots aligned over the field lines is likely to be similar 
to that of extended heating. 

Persoon et al. [1983] have reported basic transitions at 
geocentric altitudes of 1.16 /? E and 1.55 /? E in the radial 
distribution of the plasma density in the polar cap region. The 
lower-altitude transition is attributed to the transition from 
collisional to noncollisional flow. The transport process leading 
to the higher-altitude transition is not well understood. It is 
quite likely that low-level transverse heating extending in both 
time and space may be a mechanism for such a transition. We 
point out that the formation of cavities by ion heating may 
facilitate the outflow of the heavy ionospheric ions due to the 
enhanced parallel electric fields near the lower regions of the 
cavities. This is especially true if the heating and the cavity 
formation extend down to altitudes of <3000 km [Singh et al., 
1989]. 

It must be pointed out that there are numerous observations 
of transversely heated ions (see, for example, review by Klumpar 
[1986]). However, there is almost a total lack of any systematic 
observational study of how the process of transverse ion heating 
modifies the ambient plasma distribution in the magnetosphere. 
In the literature it is occasionally mentioned that ion conics are 
somehow correlated with severe plasma depletions [Klumpar, 
1979, 1981]. In this paper it is shown that plasma depletions 
are one of several consequences of transverse ion heating 
whether the heating is localized or extended. Other 
consequences, such as the density enhancement, generation of 
localized parallel electric fields, and parallel heating in front of 
the transversely heated ions, remain to be examined 
observationally. 
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Fig. 1. Geometry of the calculation. The lower and upper boundaries 
of the calculation region are at geocentric distances r=r 0 = 9371 km 
and r=r maJ =93,371 km. The distance from the lower boundary is s. 


Fig. 2. Plasma perturbations in response to a relatively weak heating 
rate of 0.24 eV/s: (a) density, (b) flow velocity, (c) parallel 

temperature, and ( d) perpendicular temperature. The curves shown 
in these panels are at times f„ = 0.922 x n for n = 1, 2, 4, and 6. 
Electron temperature T t is 1 eV; healing time Ar h is 1.844 min and 
heating region size Ar h is 210 km. 


Fig. 3. Same as Figure 2, but for a relatively strong heating rate of 2.4 
eV/s. 


Fig. 4. Electric field distribution at (a) t = /,, (A) t = t 2 , (c) t = r 4 , 
and (d) t = The solid curves are the electric field calculated solely 
from the gradients in the density profile while the circles on the curves 
are the sum of the contributions from di/ds, and dV/dt in 

equation (6). Note that the electric field is quite accurately described 
by the density gradients alone. Parameters for this figure are the same 
as in Figure 3. 


Fig. 5. Comparison between the electric fields from Figure 4 with the 
mirror force F m . The dashed curves show FJe. 


Fig. 6. Same as Figure 3, except for A/ h = 0.1844 min and heating 
rate of 24 eV/s. 


Fig. 7. Comparison between plasma perturbations for tj = 0.01 
(dashed curve) and ij = 1 (solid curve). Left- and right-hand panels 
are for early (i = tj) and late (t = t 6 ) times, respectively. Other 
parameters are the same as in Figure 3. 


Fig. 8. Same as Figure 3, but for electron temperature linearly 
increasing with altitude from 0.3 eV at s = 0 to 3 eV at s = 12,000 
km. 


Ftg. 9. Electric fields at (a) t - r lt (b) t = (c) t = r 4 , and (d) t = 

t 6 associated with the perturbation in Figure 8. 


Fig. 10. Plasma perturbations in response to impulsive heating in both 
lime and space (Ar h = 5 s, A s h = 60 km). Electron temperature T c is 
10 eV. (a, b, c, and d) Evolution of n( s), F(s), r ( (s), and 7\(s), 
respectively, (e, f, g, and h) Evolution of the electric field distribution. 


Fig. 11. Plasma perturbations in response to an extended heating 
event: (a) density, (b) velocity, (c) parallel temperature, and (d) 
perpendicular temperature. Heating rale is 0.24 eV/s. The profiles 
shown are at times r n = n x 1.844 min, n = 1, 3, 5, 7, and 9. Heating 
terminates at r = f u = 20 min. 


Fig. 12. Comparison of flow velocities (dashed curves) obtained by 
' integrating equation (10) with those (solid curves) given in Figure 116 
at times (a) r 3 , (6) r 5 , (c) r 7 , and (d) Note the good agreement 
between the two sets of curves inside the cavity. The top of the cavity 
is indicated by the arrows. 


Fig. 13. Same as Figure lid, but on a linear scale. 


Fig. 14. Pitch angle distribution at t = t 9 for the extended heating 
shown in Figure 11. 


Fig. 15. Polar wind recovery after the heating ceases at t = r n : (a) 
density, ( b ) velocity, (c) parallel temperature, and ( d) perpendicular 
temperature. 


Fig. 1. Geometry of the calculation. The lower and upper boundaries of the calculation region are at geocentric distances 
r = r Q = 9371 km and r = = 93,371 km. The distance from the lower boundary is j. 


Fig. 2. Plasma perturbations in response to a relatively weak heating rate of 0.24 eV/s: (a) density, (b) flow velocity, (c) 
parallel temperature, and ( d ) perpendicular temperature. The curves shown in these panels are at times r n = 0.922 x n for 
n = 1, 2, 4, and 6. Electron temperature T c is 1 eV; heating time Ar b is 1.844 min and heating region size Ar b is 210 km. 


Fig. 3. Same as Figure 2, but for a relatively strong heating rate of 2.4 eV/s. 


Fig. 4. Electric field distribution at (a) t = r,, (b) t = r 2 , (c) t = t 4 , and (d) t = r 6 . The solid curves are the electric field 
calculated solely from the gradients in the density profile while the circles on the curves are the sum of the contributions 
from dt/clf, dV/dr, and dV/dt in equation (6). Note that the electric field is quite accurately described by the density 
gradients alone. Parameters for this figure are the same as in Figure 3. 


Fig. 5. Comparison between the electric fields from Figure 4 with the mirror force F a . The dashed curves show FJe. 


Fig. 6. Same as Figure 3, except for Ar h = 0.1844 min and heating rate of 24 eV/s. 


Fig. 7. Comparison between plasma perturbations for r? = 0.01 (dashed curve) and tj = 1 (solid curve). Left- and 
right-hand panels are for early (r = rj and late (r = r 6 ) times, respectively. Other parameters are the same as in Figure 
3. 


Fig. 8. Same as Figure 3, but for electron temperature linearly increasing with altitude from 0.3 eV at s = 0 to 3 eV at 
s = 12,000 km. 


Fig. 9. Electric fields at (a) t = r,, ( b ) t = r 2 , (c) r = r 4 , and (d) t = i 6 associated with the perturbation in Figure 8. 


Fig. 10. Plasma perturbations in response to impulsive heating in both time and space (Ar b = 5 s, Ar b = 60 km). Electron 
temperature T e is 10 eV. (a, b, c, and d) Evolution of n(s), F(s), T^s), and TJs), respectively. ( e,f , g, and h) Evolution 
of the electric field distribution. 


Fig. 11. Plasma perturbations in response to an extended heating event: (o) density, (f>) velocity, (c) parallel temperature, 
and ( d) perpendicular temperature. Heating rate is 0.24 eV/s. The profiles shown are at times r n = n x 1.844 min, n = 
1, 3, 5, 7, and 9. Heating terminates at t = r n = 20 min. 


Fig. 12. Comparison of flow velocities (dashed curves) obtained by integrating equation (10) with those (solid curves) given 
in Figure 11 b at times (a) r 3 , (b) r 5 , (c) r 7 , and (d) Note the good agreement between the two sets of curves inside the 


cavity. The lop of the cavity is indicated by the arrows. 


Fig. 13. Same as Figure 11 d, but on a linear scale. 


Fig. 14. Pitch angle distribution at t = r, for the extended heating shown in Figure 11. 


Fig. 15. Polar wind recovery after the healing ceases at / = r tl : (a) density, (b) velocity, (c) parallel temperature, and (d) 
perpendicular temperature. 
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